Spin waves in a Bose-condensed atomic spin chain 
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The spin dynamics of atomic Bose-Einstein condensates confined in a one-dimensional optical 
lattice is studied. The condensates at each lattice site behave like spin magnets that can interact 
with each other through both the light-induced dipole-dipole interaction and the static magnetic 
dipole-dipole interaction. We show how these site-to-site dipolar interactions can distort the ground 
state spin orientations and lead to the excitation of spin waves. The dispersion relation of the spin 
waves is studied and possible detection schemes are proposed. 
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Rapid experimental progress in the realization of 
trapped degenerate quantum atomic gases has gener- 
ated fascinating opportunities to study a wide range of 
physical phenomena in atomic physics, condensed mat- 
ter physics and quantum optics. In particular, the re- 
cent success in all-optical confinement [Q, ^ and forma- 
tion H of atomic Bose-Einstein condensates provides a 
unique tool to explore the magnetic properties and spin- 
dependent dynamics of ultracold atomic gases 0] . In this 
letter, we propose a scheme to study the excitation and 
propagation of spin waves in an array of atomic spinor 
Bose-Einstein condensates confined or created in an op- 
tical lattice. 

Spin-wave phenomena play an important role in solid 
state physics ||. In solids, spin- wave excitations result 
from the exchange interaction of electrons between atoms 
in the crystal. They are usually associated with spin-1/2 
Fermi systems with an effective interaction range of a 
few angstroms, a typical lattice period in solid materials. 
Although the Bose-condensed atoms in optical lattices 
exhibit a number of close analogies to atoms in a real 
crystal, a number of differences exist. (1) The atomic 
spacing in an optical lattice is of the order of half an op- 
tical wavelength and as such is much larger than a crys- 
tal lattice period. (2) As a result, the electron exchange 
interaction is completely negligible: spin waves, if they 
exist, must be caused by other forms of long range inter- 
actions. (3) There is a large number of (bosonic) atoms 
at each lattice site, typically of the order of 1000 or more, 
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and they are subject to Bose-enhancement effects. As a 
result of these differences, the Bose condensates in an op- 
tical lattice offer a totally new environment to study spin 
dynamics in periodic structures. 

The schematic diagram of Fig. |l| shows the system dis- 
cussed in this letter. We consider for concreteness a 
one-dimensional (ID) optical lattice formed by two 7T- 
polarized laser beams counter-propagating along the y- 
axis. We assume that, in the x-z plane, Bose-condensed 
alkali atoms in their hyperfine ground state manifold are 
tightly confined by an optical dipole potential arising 
from either the transverse profile of the lattice field or 
from a separate laser. Hence a ID coherent atomic chain 
is formed along the y-axis. We employ a spinor atomic 
field theory to describe the interaction of the atoms with 
the lattice laser beams. For large detunings A — lot, — u) a 
between the frequency u>l of the laser fields and the 
atomic transition frequency io a it is possible to adiabat- 
ically eliminate the excited atomic state field operator. 
The resulting Heisenberg equations of motion for the hy- 
perfine ground-state atomic field operators ip m (r,t), in- 
cluding interatomic collisions and static magnetic dipole- 
dipole interaction, take the form || m, Bl BL 
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FIG. 1: Schematic diagram shows that a spinor condensate 
is confined in an optical lattice in the y-axis with an external 
magnetic field along the z direction. 
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where Vl(t) = U exp[-r^/W£] cos 2 (/clj/) is the light- 
induced lattice potential, with r± = \/x 2 + z 2 . The po- 
tential depth is defined as Uq = ?i|f2| 2 /6A with Q, being 
the Rabi frequency, fc^ = 2it/\l the wave number of 
the lattice beams, and Wl the beam width. The index 
to = — F, ...F denotes the Zeeman sublevels of the elec- 
tronic ground state of the atoms with angular momentum 
F. In this letter we take F = 1 for the ground state alkali 
atoms. 

The first nonlinear term in Eq. ([!]), originates from 
the photon-exchange interaction between the condensed 



atoms. It describes the light-induced dipole-dipole inter- 



action and is characterized by the quantity 
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where 7 is the single-atom spontaneous emission rate and 



e +1 )/6h denotes the dipole mo- 
ments induced by the 7r-polarized light fields. Here, i%n 
are the matrix elements of the "+" and "— " components 
of the total angular momentum operator F, and e-tio 
are unit vectors in the spherical harmonic basis. The ten- 
sor W, describing the spatial profile of the light-induced 
dipole-dipole interaction, has the form 



For F = 1, the individual condensates consist of atoms 
with three Zeeman sub-levels, hence they behave as col- 
lective spin magnets in the presence of external magnetic 
fields or spin-dependent interactions. Such spin magnets, 
localized along the lattice axis, form a ID coherent Bose 
atomic spin chain. Using Eq. ([!]) and Eq. ([!]), and ignor- 
ing both the non-resonant and spin-independent constant 
terms, we can construct the Hamiltonian describing this 
spin chain, 
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where 11 is the unit tensor, r = r/|r| and £ — /ci|r|. 

The two-body ground-state collisions and magnetic 
dipole-dipole interaction are described by the potentials 
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respectively, where A s and A a are related to the s-wave 
scattering lengths of the spinor condensate (Tcfl. 

The optical potential associated with a sufficiently 
deep optical lattice is equivalent to a periodic array of in- 
dependent "microtraps" ® . If the depth of each of those 
is large enough, an array of independent Bose conden- 
sates can be formed in the lattice, and it is convenient to 
expand the spinor atomic field operator as 



also introduced an external magnetic field B = B^eo 
whose strength is strong enough to polarize the ground 
state spin orientations of the atomic chain along the 
quantization axis zQ]. The parameter 75 = —/j,b9f is 
the gyromagnetic ratio, /is being the Bohr magneton and 
gp the Lande (/-factor. 

The first term in Hamiltonian (||) results from the spin- 
dependent interatomic collisions at a given site, with 
\' a = (l/2)A a /d 3 r|4(r)| 4 . The last two terms de- 
scribe the site-to-site spin coupling induced by both the 
static magnetic field and light-induced dipolar interac- 
tions. The coupling coefficients have the explicit forms 



]fa(r)a m (i) 
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where fa is the condensate wave function for the ith 
microtrap and the operators a m (i) satisfy the bosonic 
commutation relations [a m (i), a„(j)] = S mn 6ij. For deep 
enough microtraps the spatial overlap between the indi- 
vidual condensate wave functions is negligible, and they 
can be considered as independent. Under this tight- 
binding condition, the spatial wave function of the ith 
condensate is then determined by the Gross-Pitaevskii 
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equation 
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fiifa(r), with Vi{v) being the potential near the z-th mi- 
crotrap and Ni = ^ m {d m {i)a m {i)) the number of con- 
densed atoms at the site, and Hi the chemical potential. 
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where we have introduced a cut-off function f c (r) — 
cxp(— r 2 /L 2 ) to describe the effective interaction range 
of the light-induced dipole-dipole interaction, L c = _/V 7 /c 
being the coherence length associated with the collective 
spontaneous emission of N atoms |i~H . 

The physics implicit in Hamiltonian (j^) is quite clear. 
Consider first the situation without site-to-site coupling, 
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Jij — 0. In the presence of a sufficiently strong exter- 
nal magnetic field, the spins align themselves along the 
quantization axis z. For ferromagnetic condensates as in 
the case of 87 Rb (A' Q < 0)|Q — or in the presence 
of a strong external magnetic field for anti-ferromagnetic 
condensates — the ground state of the Hamiltonian is 
\GS) = \N,-N), where N = Y<i N i is thc total atomic 
number in the lattice. The total spin at site i has the 
expectation value (5f ) = —Nih, where the factor Ni is 
due to Bose enhancement. 

For ^ 0, the situation changes drastically: thc 
transfer of transverse spin excitation from site to site is 
allowed, resulting in the distortion of the ground state 
spin structure. This distortion can propagate from site 
to site and hence generate spin waves along the Bose con- 
densed atomic spin chain. From Hamiltonian (^), we can 
derive the Heisenberg equations of motion for the spin 
excitations as 
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where we have replaced the spin operator S% by its 
ground state expectation value in the mean-field ap- 
proximation. The frequencies luq = —"/bB and Au) q = 
^^Zj^qJqjNjh describe the precessing of the gth spin 
caused by the external magnetic field and the static mag- 
netic dipolar interaction. The site-to-site spin coupling 
coefficients \qj = 2J q jN q h determine the propagation 
of the spin waves. From Eq. (||), we observe that the 
light-induced dipolar interaction only contributes to spin 
coupling in the x-y plane. This is because the 7r-polarized 
lattice beams only induce an effective dipole moment in 
that plane. 

To further understand how Eq. (Q) determines the exis- 
tence and propagation of spin waves, it is helpful to con- 
sider the special case where the lattice is infinitely long 
and thc spin excitations arc in the long- wavelength limit. 
We can then reexpress Eq. (R) in its continuous limit by 



the replacements S, 
(uj + Auj q ) -> u)(y), 
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s(y,t), Xq 3 -> x(y - y'), and 



l ^JT^ = u{v)s{y > t] 'Y L l dy ' xiv ~ y ' )S{y ' 1 1] ' (8) 

In the long- wavelength limit, the last integral term in 
Eq. (|^) can be evaluated up to the second-order expan- 
sion of S(y',t). This leads to an effective Schrodinger 
equation 



. dS(y,t) 
'' dt 



' 2 dy 2 



S(y,t) 
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where we have defined (5 n = (2/Al) / drjx{ r l)' r l 2n f° r 
n = 0,1. Evidently, from Eq. (0), S(y,t) describes the 
"waves" caused by spin excitations in the x-y plane. 
These waves can be associated with the center-of-mass 
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FIG. 2: Magnon dispersion spectrum. In the calculation, we 
have assumed that the spatial dimension along y-axis of the 
condensate in each lattice site is much less than Al, while 
in the transverse (x-z) plane, the condensate has a Gaus- 
sian shape with a width w. We have taken \l = 1/xm and 
7|f2| 2 /A 2 = 10 3 . We have used a total number of 100 lattice 
sites with 2000 atoms in each site. 



motion of a quantum-mechanical particle of effective 
mass m = Being similar to the phonon associ- 

ated with sound waves, the excitations associated with 
the spin waves are usually referred to as "magnons" |g . 

We can evaluate the magnon dispersion rela- 
tion by solving the discrete wave equation (Q) for 
waves of the form S q ^ = dk exp[— i(u>o + Aluq — 
2 Sj>o^Oj)i]exp[i(fc(7A i /2 - u k t)}. For simplicity, we 
assume a lattice with each site having the same number 
of atoms, Ni = Nq. This yields the magnon dispersion 
relation 



= 2^xoi[l-cos(jfcA i /2)] ! 



(10) 



j>0 



where we have taken q = in the coefficients of Eq.(|7]) 
since their values are independent of q in the case at hand 
as long as we have a sufficient number of lattice sites. 

Figure || shows the magnon dispersion spectrum and 
its dependence on the transverse width w of the Bose 
condensates. As w increases the dipolar interactions 
among atoms tend to cancel each other, thereby reduc- 
ing the excitation frequency of the spin waves. As a re- 
sult, it becomes difficult to excite the spin waves in the 
coherent atomic chain through the dipole-dipole interac- 
tion for condensate widths much larger than one optical 
wavelength. In addition, our numerical calculations show 
that for typical experimental parameters the strength of 
the light-induced dipolar interaction dominates over the 
magnetic dipolar interaction in the determination of the 
coupling coefficients \0j- 

Having established the existence of spin waves in con- 
densate lattices, the question remains to determine how 
to detect them. Any optical or magnetic method which 
can excite the internal transitions between the atomic 
Zeeman sublevels can be used for this purpose. On the 
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m=- 1 

FIG. 3: Magnon detection schemes. 



other hand, the detection scheme should be able to distin- 
guish different spin wave modes (see discussion below) . A 
natural choice consists in employing Raman transitions, 
as shown in Fig. ||. In this case the combination of a 
circularly-polarized and a 7r-polarized Raman beam will 
couple all three Zeeman sublevels. This exactly creates 
the spin transition associated with the operators 
The existence of spin waves can then be detected by mea- 
suring the absorption of one of the Raman beams. The 
absorption spectrum is proportional to the sum of the 
transition probabilities among the Zeeman sublevels, 



P (X 
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| ($ k \S+\GS)\ 2 sin 2 [(^-^-^)r/2] 
2 [{v - uj p - 0J k )r/2} 2 



where the kets |$fe) denote excited states with spin waves 
of frequency lou, v is the frequency difference between 
the two Raman beams, r is the measurement time and 
ujp = ujo+Aujq — 2 J2j >Q Xoj defines the total spin process- 
ing frequency. Absorption resonances occur whenever the 
frequency v is tuned to match a spin- wave frequency u>k- 
For an infinitely long lattice, this will produce a broad- 
ened absorption spectrum whose width characterizes the 
existence of spin waves. In practice, though, the optical 
lattice is finite with length L. Such lattices only allow the 
excitation of spin waves with discrete wave numbers re- 
sulting from the resonance excitation condition kL = mr, 
k being the wave number of the spin waves in Eq. (10). 
As a result, the absorption spectrum will exhibit a multi- 
peak structure. In the long-wavelength limit, the interval 
between peaks is proportional to 



L j>o 



with Nl = L/(\l/2) being the total number of lattice 
sites. 

In current experiments in optical lattices, Nl is in the 
range of 10 ~ 100, and each lattice site can accommo- 
date a few thousand atoms, in which case X^>o Xojj 



can reach values of about a few MHz. This leads to a 
requirement for the frequency measurement precision of 
about 10 ~ lOOKHz. This is achievable with current 
techniques. 

Alternatively, one can also carry out measurement in 
momentum space. The magnons associated with spin 
waves of wave number k have momenta p = hke y . If 
one observes the Raman scattering using two Raman 
beams through the Bose gas, the momentum conserva- 
tion between the magnons and Raman photons requires 
Ak = ke y with Ak being the difference of wave vectors 
between two Raman beams. Hence the momentum dis- 
tribution of the scattered Raman photons can identify 
the existence of the spin waves. 

In conclusion, Bose condensates in an optical lattice 
offer a new tool and test ground to study the quantum 
spin phenomena. The Bose statistics in each lattice site 
makes the atoms coherently behave as a spin magnet. 
Such an array of spin magnets not only exhibit fascinat- 
ing spin dynamics as shown in this letter, but also may 
find potential applications in quantum information and 
computation. This opens new opportunities for future 
research. 
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